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The Forward-Backward Asymmetry (AFB) in Z' physics is commonly only perceived 
as the observable which possibly allows one to interpret a Z' signal by distinguish¬ 
ing different models of such (heavy) spin-1 bosons. In this paper, we examine the 
potential of AFB in setting bounds on or even discovering a Z' at the Large Hadron 
Collider (LHC) and show that it might be a powerful tool for this purpose. We anal¬ 
yse two different scenarios: Z'-bosons with a narrow and wide width, respectively. 

We find that, in the first case, the significance of the AFB search can be comparable 
with that of the bump search usually adopted by the experimental collaborations; 
however, being a ratio of (differential) cross sections the AFB has the advantage of 
reducing systematical errors. In the second case, the AFB search can win over the 
bump search in terms of event shape, as the structure of the AFB distribution as 
a function of the invariant mass of the reconstructed Z'-boson could nail down the 
new broad resonance much better than the event counting strategy usually adopted 
in such cases. 


I. INTRODUCTION 

Heavy neutral Z'-bosons arise in a number of theories that extend the Standard Model (SM) 
gauge group by adding an extra U{1) symmetry. The most common Z'-boson benchmark 
models can be divided in three main classes: Eq models, Generalized Left-Right (GLR) 
symmetric models and Generalized Standard Models (GSM), see, e.g., the review in [1] and 
references therein. All these models predict a relatively narrow width for the Z'-bosons, so 
that Vz'/Mz' varies in the 0.5 — 12% range. The lowest Vz' j^z' value is realised in the 
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model from the class while the biggest value appears in the Q-model belonging to the 
GSM class. 

Experimental searches for a heavy Z'-boson at the LHC are usually interpreted in the context 
of the Sequential Standard Model (SSM), which is part of the GSM class [2]. This benchmark 
scenario just includes one extra neutral vector boson with couplings to fermions identical 
to those of the corresponding SM Z-boson and no mixing with the neutral Electro-Weak 
(EW) SM boson. Being nothing but a heavier copy of the SM Z-boson, this Z'-boson is 
characterized by a narrow width: Vz' jMz' — 2.8%, including the Z'-decay into top-antitop 
pairs above threshold. Dedicated search strategies therefore assumes that the new heavy 
resonance is narrow and can be described by a Breit-Wigner line-shape, standing over the SM 
background, when looking at the distribution in the invariant mass of the Z'-boson decay 
products. In this way, the new physics signal is thought to have a well dehned peaking 
structure, concentrated in a small interval centred around its mass. On the basis of this 
assumption, the 95% Gonhdence Level (G.L.) upper bound on the cross section is derived 
and limits on the mass of the Z'-boson are extracted within the above mentioned benchmark 
models. In the case of a narrow width Z' scenario, even interference effects can be accounted 
for without substantially altering the described experimental approach [3] . 

However, the narrow width hypothesis is quite strong, even if well motivated. There 
exist in fact counter examples of theories where the predicted Z'-boson is characterized by 
a large width. For example, this can be realized in Technicolor |3] scenarios, Gomposite 
Higgs Models jS] or in more generic setups where the Z'-boson couples differently to the 
hrst two fermion generations with respect to the third one PEI or else interacts with the 
SM gauge bosons in presence of mixing [8], so that large Vz'lMz' values are induced by 
the additional Z' decay channels which then onset in all such cases. As a consequence, the 
resulting resonance is wide and, instead of having a well dehned Breit-Wigner line shape, 
it will appear as a broad shoulder spreading over the SM background. In such models, 
the ratio T^z/M^/ can easily reach the 50% value or higher making a Breit-Wigner line 
shape based analysis dehnitely inappropriate. The experimental collaborations have not 
tackled this kind of scenarios yet and respective dedicated search strategies for wide Z' 
particles are absent at the moment. The only frameworks which have been analyzed in 
order to interpret possible non-resonant deviations from the SM are the graviton production 
within the ADD model P EO] and the contact interaction within the left-left isoscalar model 
mm- For details on the parametrization of the differential cross section within these two 
scenarios we refer to Ref. [I3] and references therein. Both processes might give rise to an 
excess of events spread over the SM background. In this ‘effectively’ non-resonant case, the 
experimental analyses are essentially counting experiments: an excess of events searched out 
of an estimated SM background. To make the analyses more robust, the same background 
is often estimated with multiple data-driven methods. Kinematical cuts are then optimized 
in order to maximize the discovery/exclusion potential at the LHG. Despite this, as one can 
understand, this analysis can be fragile. The experimental results heavily rely on the good 
understanding and control of the SM background, as the new physics signal is not expected 
to have a dehnite shape. The choice of the control region, needed to dehne the functional 
form of the SM background to be used in the regions where there might be some signal, is 
not trivial. Interference effects between new physics signal and SM background can indeed 
affect the low scale region of the distribution in the invariant mass of the Z'-boson decay 
products, proving the assumption that the control region is new physics free to be simply 
false. Under these premises the experimental analyses could get quite complicated. On the 
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one hand, the possible presence of wide objects could in fact shrink the new physics free 
region, owing to an increase of the interference effects driven by the large width. On the 
other hand, the wide Z'-bosons could easily escape detection in the bump searches due to 
the same interference effects combined with the absence of a resonant peaking structure in 
the di-lepton invariant mass distribution. All this then conspires to make the discovery of a 
wide Z' very problematic. 

In this paper, we study the potential of a complementary observable in order to probe new 
massive Z' objects either narrow or broad. The variable we analyse is the Forward-Backward 
Asymmetry (AFB). In the literature, this observable is usually advocated in the second stage 
of the data analysis process to interpret experimental results after a possible discovery of a 
new spin-1 particle using a standard bump hunt. The AFB is indeed used to disentangle 
different models with Z'-bosons and nail down the underlying theory. This procedure relies 
on the assumption that the new heavy Z'-boson is characterized by a narrow width and it 
would be discovered via the bump search. Our purpose is to show that the AFB can be used 
not only for interpreting a possible discovery but also in the very same search process. We 
show that the AFB observable can be associated to the default resonance search to improve 
and/or extend the discovery potential for both narrow and wide Z's. Focussing on AFB, 
we aim at establishing the methods needed to study Z'-boson production at the CERN 
LHC in the Drell-Yan channel, giving rise to di-lepton pairs in the final state: pp —)■ 
with / = e, /i. This production process is particularly clean and thus represents the golden 
channel for Z'-boson discovery at the LHC. 

The paper is organized as follows. In Section]^ we derive current and projected bounds 
for Z' model benchmarks for the LHC at 8 and 13 TeV, respectively, for the models presented 
in Ref. [1]. In Section HI, we discuss the role of the forward-backward asymmetry within 
the Z’ physics, its reconstruction and statistical uncertainty. We also discuss the effect of 
a rapidity cut on the di-lepton system upon the signal signihcance. This cut is commonly 
implemented in order to increase the efficiency in guessing the quark direction in pp collisions, 
which is needed to reconstruct the AFB observable. The drawback of applying it is a decrease 
of the number of signal events with the consequent depletion of the AFB signihcance. The 
outcome is that this stringent cut can be relaxed for the range of Z'-boson masses which 
we will be looking for in the next LHC run. In Section |IV[ we discuss the role of the AFB 


in searches for narrow width Z'-bosons and systematically analyse Z' model benchmarks 
confronting AFB and bump search. We will show that the signihcance from the AFB can 
be comparable with that obtained from the cross section studies over the same invariant 
mass distribution of the di-lepton system. In this scenario, the advantage of using the AFB 
observable would consist in minimizing the systematics, as the AFB is a ratio of (diherential) 
cross sections. Moreover, the two observables (cross section and asymmetry) could be of 
mutual support to make the claim of a possible new physics discovery more robust, if the 
bump search itself would provide only a mild evidence for a Z' state. 

In Section |V| we analyse the role of AFB in searches for wide Z' particles. We consider 
two benchmark models which predict a wide Z'-boson with ratio Tz'/Mz' of the order of 
several tens of percent. In this case, again, the AFB can be complementary to the resonance 
search and have a distinctive line shape contrary to the invariant mass distribution of the 
cross section, which could well mimic the background shape. That is, the latter would only 
give rise to an excess of events evenly spread over the SM background, which is of difficult 
interpretation and in many cases of difficult measurement owing to uncertainties in the 
background modelling. Finally, in Section [VH| we summarize and conclude. 
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FIG. 1: (Colour online) (a) 95% C.L. upper bound on the Z'-boson production cross section 
in Drell-Yan normalized to the SM cross section on the Z-boson peak: Rfj = cr{pp —)• Z' —)• 
/ (t{pp — 7- Z, 7 — )• /■*“/“) with I = e, p. The combined analysis of the di-muon and di-electron 
channels has been produced by the CMS collaboration with a data sample collected at the 8 TeV 
LHC, corresponding to an integrated luminosity of 20.6 and 19.7 fb~^ respectively [13]. Theoretical 
predictions for the class of the Eq models are superimposed to extract the corresponding Z'-boson 
mass limits. As described in the text, in order to match theoretical predictions and experimental 
results, the optimal cut on the invariant mass of the di-lepton pairs has been implemented: AM = 
\Mii — Mz'\ < 0.05 Elhc for Elhc = 8 TeV. |(b)| Same for the other two classes of GSM and 
GLR models. 


II. BOUNDS ON THE Z’-BOSON MASS 

In this section, we re-obtain the existing bounds on the mass of the Z'-boson from the 
aforementioned benchmarks models as obtained by, e.g., the CMS collaboration, after the 7, 
8 TeV run with about 20 fb~^ of accumulated luminosity, assuming a narrow width: these 
can be found in Ref. jTHj. After conhrming current CMS limits obtained under the above 
assumption, we further present those obtained by taking into account the full width effect 
as well the intereference corrections. Finally, we produce projection limits for LHC Run 2. 

We therefore start by scanning over the thirteen benchmark models predicting a Z'-boson 
characterized by a narrow width (Tz'/Mz' < 5%) summarized in Ref. [1], and extract the 
limits on by making use of the 95% C.L. upper bound on the Z'-boson production cross 
section in Drell-Yan, a{pp —)■ Z' —)■ e~^e~, p~). In order to reduce systematic uncertainties, 
the experimental analysis normalizes the Z'-boson production cross section in Drell-Yan to 
the SM Z-boson cross section on peak. As shown in Fig. 1, the 95% C.L. upper bound is 
indeed given on the ratio Ro- = a{pp —)■ Z' —)■ e^e~, p~)/a{pp —>■ Z, 7 —>■ e+e“,/r+/r“). 
The use of this ratio in fact cancels the uncertainty in the integrated luminosity and 
reduces the dependence on the experimental acceptance and trigger efficiency. 

We calculate this ratio, R^, at the Next-to-Next-to-Leading Order (NNLO) in QCD using 
the WZPROD program [T411I6] (which we have adapted for Z' models and new Parton 
Distribution Function (PDF) sets [I]) and the CTEQ6.6 package [T7] . 

The NNLO QCD contributions give rise to a R-factor which depends on the energy 
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2700 

2560 

2620 

2640 

2600 

2570 

3040 

2950 

2765 

3260 

2900 

3135 

3720 


TABLE I: Bounds on the Z'-boson mass derived from the latest direct searches performed by 
CMS at the 8 TeV LHC with integrated luminosited L = 20/6“^. We consider thirteen different 
models with an extra U'{1) gauge group predicting a new heavy neutral boson characterized by a 
narrow width. From left to right, the columns indicate the Mz' limit in GeV within the Eq, GLR 
and GSM class of models. 

scale, thus we fully take into account such a dependence. The NNLO prediction for the 
SM Z,'j production cross section, a{pp —)• Z, 7 —/+/“) with I = e,fi, in the mass window 
of 60 to 120 GeV is 1.117 nb. With all these ingredients at hand, we compute as a 
function of the mass of the new heavy Z'-boson, and derive the corresponding limits 

for all benchmark models. Fig. la shows the bounds on all Eq models, while Fig. lb 
displays the results for the remaining two classes of models, GLR and GSM. As previously 
mentioned, traditional experimental analyses work under the hypothesis that the signal has 
a Breit-Wigner line shape and performs the analysis in a restricted search window around 
the hypothetical mass of the Z'-boson. This approach is theoretically motivated by the 
benchmark models, all predicting a narrow width Z'-boson, and by the will to perform an 
analysis as much as possible model independent. One should stress that the GMS analysis 
makes use of a dedicated cut on the invariant mass of the di-lepton pairs: \Mii — Mz'] < 
0.05 X F^lhc where F^lhc is the collider energy. This cut was designed so that the error 
in neglecting the (model-dependent) Finite Width (FW) and interference effects (between 
7 , Z, Z') are kept below 0(10%) for all models and for the full range of allowed Z' masses 
under study, thus following the recommendations of |3]. This procedure thus continues to 
allow for a straightforward interpretation of the extracted mass bounds in the context of any 
theory predicting a narrow Z'-boson. At this stage of our own analysis, we work under the 
very same setup, for validation purposes, with the notable exception that we allow for the 
aforementioned FW and interference effects, unlike the experimental results which assume 
the so-called Narrow Width Approximation (NWA), wherein the (narrow) Z' is actually 
produced on-shell. Table [T] summarizes the bounds we obtain. They reproduce the GMS 
limits very well in general, within the accuracy of 1 — 2%. The only slight exception is the 
Q-model in the GSM class where our limit, based on the present analysis accounting for full 
width and interference, is different from the GMS one by about 5%. Yet, this is well in line 
with expectations, as this model is the one yielding the largest width. 

But let us set the stage in some more detail. In Fig. 2 we show the behaviour of the 
new physics signal for two representative scenarios: the model (Fig. 2a) and the SSM 
benchmark scenario (Fig. 2b). The solid line represents the full new physics signal, that 
is, Z'-boson production and decay including the interference with the SM background. The 
dashed line gives instead the pure Z'-boson signal, neglecting the interference. (As evident 
from the plots, in both cases we allow for FW effects of the Z'-boson.) As one can see, 
the shape of the distribution in the invariant mass of the di-lepton system is quite model 
dependent off peak. This is due to the fact that it can be heavily affected by model dependent 
interference effects there. The sign of the interference is not dehned a priori. It can be either 
positive, like in the Eq models, or negative, like in the two other classes of GLR and GSM 
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(b) 


(a) Differential cross section in the invariant mass of the di-lepton system coming from the 

We consider the model. 


FIG. 2: 

Z'-boson production and decay in the Drell-Yan channel: pp e"’“e 
The solid line shows the complete new physics contribution to the invariant mass distribution, that 
is the Z' signal plus the interference with the SM background. The dashed line represents instead 
the pure Z' signal, neglecting the intereference. The plot has been produced for the 8 TeV LHC, 
no cuts are applied. |(b)| Same for the SSM benchmark scenario. 


models. In addition, its magnitude can be quite sizeable. The off peak tail of the signal 
distribution in the di-lepton invariant mass is thus highly model dependent in the low mass 
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FIG. 3: (Colour online) (a) Discovery potential of the 13 TeV LHC for the Eq class of models. 
We plot the 5a contours as a function of Z'-boson mass and luminosity. We perform a combined 
analysis over e'''e“ and pairs and assume the Axe factor given by CMS at the 8 TeV LHC. 

|(b)| Same for the GSM and GLR classes of models. 


region, leading to either an excess or a depletion of the total number of expected events as 
compared to the SM background, according to the sign of the interference. 

Furthermore, it is clear that the fully integrated cross section for the complete new 
physics signal, that is, Z'-boson production and decay including the interference with the 
SM background, is not always a uniquely dehned variable. In the SSM, and more generally 
in all GSM and GLR models, the signal can manifest itself as a negative correction to the 
differential cross section (solid line) at low masses. Similarly, the fully integrated cross 
section for the pure Z'-boson production and decay, neglecting the interefences, can also 
give quite a wrong picture. Taking into account the low mass tail of the invariant mass 
distribution can overestimate the Z'-boson signal by a large factor. For a SSM Z' with 
mass Mz' = 3 TeV, we have that the fully integrated cross section for the pure signal is 
az' = 0.17/fe while the complete signal cross section, integrated in the mass window where 
it is positive dehnite, is equal to az' + Interference = 0.06 fb. In this case, taking into 
account the unphysical tail (which in the SSM is in reality washed out by the destructive 
interference) leads to overestimate the Z' signal by a factor of 3 and consequently to extract 
more stringent limits that are erroneous. In essence, a shape analysis of the signal over the 
full invariant mass region is very challenging. Thus, the dehnition of the observable to be 
used to interpret the data and extract the mass bounds on the hypothetical Z'-boson must 
indeed be appropriately chosen. 

However, thanks to the approach recommended in j3], all such extreme effects are avoided 
and, in the instance, we can conclude that our code for the simulation of Drell-Yan processes 
which might receive a contribution from a narrow Z'-boson exchange, pp —?• 7 , Z, Z' —)■ 
has been validated against the GMS results. In short, by hnally taking into 
account the published acceptance x efficiency corrections {Axe), we can indeed reproduce 
the above mentioned extracted bounds on the Z'-boson mass by self-consistently evaluating 
signal and background. In doing so, we have applied the Poisson statistics for computing 
the signihcance and hence the exclusion limits on Mz'- 

By using such a code, in Figs. 3 and 4, we now project discovery and exclusion potential 
of the upgraded LHG, which will run at 13 TeV. The hrst two plots show the LHG discovery 
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TABLE II: Projection of discovery limits (first row) and exclusion limits (second row) on the Z'- 
boson mass from direct searches at the forthcoming Run II of the LHC at 13 TeV. We assume the 
original design value for the integrated luminosited: L = 300/6“^. We consider thirteen different 
models with an extra U'{1) gauge group predicting a new heavy neutral boson characterized by a 
narrow width. From left to right, the columns indicate the Mz' limit in GeV within the Eq, GLR 
and GSM class of models. 

potential for the Eq models (Fig. 3a) and for the remaining other two classes GLR and GSM 
(Fig. 3b). Figs. 4a and 4b display instead the LHG exclusion potential for Eq, GSM and 
GLR models, respectively. In deriving these results, we stress that we have included the Axe 
factor extracted by the analyses performed by GMS at the 8 TeV LHG, thereby implicitly 
assuming that no signihcant departures in this respect are assumed at the upgraded GERN 
machine. These projections are valid only for narrow width Z'-bosons for which the optimal 
observable is the invariant mass of the di-lepton system, used in the standard bump search 
performed by the experimental collaborations. From the above plots, one can conclude that 
the 13 TeV run of the LHG should be able to discover a Z'-boson with mass up to about 
4500 GeV and 5600 GeV within the Eq and GSM/GLR class of models, respectively. If 
nothing is found, the exclusion limits will be pushed up to 5300 GeV and 6400 GeV for an 
Eq Z'-boson and a GSM/GLR Z'-boson, respectively. Table [H| summarizes the LHG Run 
H potential for discovering or excluding a Z'-boson of a certain mass within all considered 
thirteen models. These projections have been obtained for the design value of the integrated 
luminosity: L = 300/6“^. 

This concludes the section on the state-of-the-art of narrow width Z'-bosons and their 
search using traditional methods. 


III. THE FORWARD-BACKWARD ASYMMETRY 


In this section, we dehne the forward-backward charge asymmetry (AFB) and we discuss 
its role in Z'-boson searches other than data interpretation. In the literature, the AFB has 
been long exploited to help disentangling the various theories predicting an extra heavy 
neutral boson and tracing back the Lagrangian parameters (see, for example, [IH] and [12] 
and references therein). This is not an easy task and the sensitivity of AFB measurements 
to new physics like additional Z'-bosons has therefore received a lot of attention in the past 
years. For Drell-Yan processes, AFB is dehned from the angular distribution 


da 

d cos 6* 


spin,col 



3 


X; + C0s-^»nc| + 2cos»;C«| 


( 1 ) 


where 0 / is the lepton angle with respect to the quark direction in the di-lepton center- 
of-mass frame (GM), which can be derived from the measured four-momenta of the di- 
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(a) 



(b) 


FIG. 4: (Colour online) (a) Exclusion potential of the 13 TeV LHC for the Eq class of models. 
We plot the 5a contours as a function of Z'-boson mass and luminosity. We perform a combined 
analysis over e'''e“ and pairs and assume the Axe factor given by CMS at the 8 TeV LHC. 

|(b)| Same for the GSM and GLR classes of models. 


lepton system in the laboratory frame. The AFB is indeed given by the coefficient of the 
contribution to the angular distribution linear in cos^^*. In Eq. V§ is the invariant mass 
of the di-lepton system, and Pi and Pj are the propagators of the gauge bosons involved in 
the process. At the tree-level, the Drell-Yan production of charged lepton pairs is mediated 
by three gauge bosons: the SM photon and Z-boson and the hypothetical Z'-boson. These 
three vector boson exchanges all participate in the matrix element squared. We thus have: 


P P 


{§ - Mf){s - Mf) + 

{{s - MfY + MfVj) {{s - Mfy + MfT]) 


( 2 ) 


where Mi and Tj are the mass and the width of the gauge bosons involved and zA = 
{ 7 , Z, Z'}. Finally, the factors Cg and in the angular distribution given in Eq. ^ are 
the parity symmetric and anti-symmetric coefficients which are functions of the chiral quark 
and lepton couplings, and to the z-boson with z = { 7 , Z, Z'}: 


Cs = (qWl + QRQR){elei + 44 ) > 
~ (44 ~ 44 ) (44 “ 44 )- 


( 3 ) 

( 4 ) 


One can conveniently compute the forward (F) and backward (B) contributions to the total 
cross section integrating over opposite halves of the angular phase space: 


dap = 


dds = 


-d cos e* = + 8&i), 


Iq d cos 6 1 
da 


1287r ^1 + 5, 




dcosO'^ 


s _ p p 

dcosor = Y 


1287r ^1 + 5, 


( 8 ) 

( 6 ) 




where z and j sum over the mediating resonances, { 7 , Z, Z'}. 

From the above expressions one can immediately see that the total cross section, a = ap+as, 
depends uniquely on the parity symmetric coefficient Cs- Conversely, the difference between 
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forward and backward cross sections, ap — ctb, preserves only the contribntion proportional 
to the parity antisymmetric coefficient Ca- This is the term which is related to the AFB. 
One can thns dehne the AFB as the difference between forward and backward cross sections 
normalized to the total cross section: 


da 


dap + das = , 


128 ^ 1 + 6 ; 


V 


ddp — das 8 s PiPj 
dap + dds 37 r(j ^ 1 + 5ij ^ 


(7) 


with the SM backgronnd corresponding to ij = 77 , ZZ, yZ and the new physics given 
by ij = 7 Z', ZZj Z'Z'. In the light of the above discnssion, total cross section and AFB 
depend on different combinations of Z'-boson conplings to ordinary matter. For that reason, 
the AFB can give complementary information abont the strnctnre of snch conplings when 
compared to the total cross section. This featnre has motivated several anthors to stndy the 
potential of the AFB observable in interpreting a possible Z'-boson discovery obtained in 
the nsnal resonance hnnt as in Ref. [T8ti20] . Onr point is that AFB can also be a powerfnl 
tool to search for new physics. 


A. The reconstructed AFB 


The AFB is obtained by integrating the lepton angnlar distribntion forward and backward 
with respect to the qnark direction. As in pp collisions the original qnark direction is not 
known, one has to extract it from the kinematics of the di-lepton system. In this analysis, 
we follow the criteria of Ref. |2I] and simnlate the qnark direction from the boost of the di¬ 
lepton system with respect to the beam axis (^-axis). This strategy is motivated by the fact 
that at the pp LHC the di-lepton events at high invariant mass come from the annihilation of 
either valence qnarks with sea antiqnarks or sea qnarks with sea antiqnarks. As the valence 
qnarks carry away, on average, a mnch larger fraction of the proton momentnm than the sea 
antiqnarks, the boost direction of the di-lepton system should give a good approximation 
of the quark direction. A leptonic forward-backward asymmetry can thus be expected with 
respect to the boost direction. In contrast, the subleading number of di-lepton events which 
originate from the annihilation of quark-antiquark pairs from the sea must be symmetric. 

As a measure of the boost, we dehne the di-lepton rapidity 


Vii 



E + P, 
E-P, 


( 8 ) 


where E and P^ are the energy and the longitudinal momentum of the di-lepton system, 
respectively. We identify the quark direction through the sign of zhj. In this way, one can 
dehne the reconstructed forward-backward asymmetry, from now on called A*p^. Namely, 
we have dehned A*pp using the 61 reconstructed angle, which is the angle between the hnal 
state lepton and the incoming quark direction in the center-of-mass of the di-lepton system. 
As the AFB reconstruction procedure relies on the correlation between the boost variable, 
i/ii, and the direction of the incoming valence quark, it is therefore more likely to pick up the 
true direction of the quark for higher values of i/if. Increasing the probability of identifying 
the direction of the quark would lead to an observed value of A*pp that is closer to the ‘true’ 
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(a) 



(b) 


FIG. 5: (Colour online) |(a)p pper plot: Probability of getting the correct direction of the valence 
up-quark at the 13 TeV LHC via the boost direction of the di-lepton system, given by the sign of 
the di-lepton rapidity yij, as a function of the modulus \yii\ = \yuu\ for six different invariant mass 
windows scanning from 500 GeV to 5000 GeV and beyond. Lower plot: Differential luminosity as 
a function of \yii\ for the correctly assigned quark pair (dashed-line) and for the full sample (solid 
line). |(b)| Same as (a) for valence down-quarks. 


value of ApB if one were able to access the partonic CM frame. The tradeoff occurs in the 
reduction of statistics which impacts the signihcances the other way. The general dehnition 
of signihcance S between predictions of an observable O with uncertainty 60 from two 
hypotheses is 

^ _ \Oi — Q2I 

“ ^ysoj + doi 

The statistical uncertainty on the AFB is given by 


SApn — 


(TfCTb 


C {ap + aB^ 




aC 


N 


where C is the integrated luminosity and N the total number of events. One can thus see 
that the signihcance is proportional to the root of the total number of events. Imposing a 
stringent cut on the boost variable, yij, would then improve the reconstructed AFB guiding 
it towards its true line shape, but it will decrease the statistics. In the next subsection, the 
subtle balance between line shape gain and statistics loss in maximizing the signihcance via 
the di-lepton rapidity cut will be discussed in detail. 


B. On the di-lepton rapidity cut 

As discussed in the previous section, since the true quark direction is not known in pp 
collisions, at the LHC one has to extract it from the kinematics of the di-lepton system. In 
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FIG. 6: (Colour online) Upper plot: Probability of getting the correct direction of the valence 
quarks at the 13 TeV LHC via the boost direction of the di-lepton system, given by the sign of the 
di-lepton rapidity Uif, as a function of the di-quark (or di-lepton) invariant mass for five different 
cuts on the di-lepton rapidity. Lower plot: Differential luminosity as a function of the di-lepton 
invariant mass for the correctly assigned quark pair (dashed-line) and for the full sample (solid 
line). Here we average on both up and down-quarks. 


this paper, the valence quark direction is approximated by the boost direction of the 
pairs with respect to the beam axis, that is given by the sign of the di-lepton rapidity Uii 
dehned in Eq. The correctness of this assignment as a function of yn has been studied 
in Ref. |2I] for di-lepton events with invariant masses above 400 GeV. In this section, we 
further analyse this issue by investigating the energy scale dependence of the probability of 
getting the true quark direction via the sign of yij. In Figs. 5a and 5b, we plot the fraction of 
events with the correctly assigned direction for up-quarks and down-quarks, respectively, as 
a function of \yii\ for different invariant mass windows of the di-lepton system. The fraction 
of correctly assigned events increases with the rapidity, conhrming the results presented in 
literature [21]. The additional information contained in Fig. 5 is that such an increase 
depends on the energy scale. For di-lepton invariant masses of TeV order the probability 
of getting the true quark direction becomes more than 80% for a rapidity cut \yij\ > 0 . 8 . 
For higher invariant masses, beyond the present Z'-boson limits of 0(3 TeV), the same 
probability can be obtained by imposing a lower rapidity cut: \yii\ > 0.35. Moreover, up- 
quarks and down-quarks respond differently to the \yii\ cut. The probability of getting the 
correct direction is higher for up-quarks than for down-quarks, at hxed \yij\ value. In Fig. 6 , 
the fraction of correctly assigned events is shown as a function of the invariant mass for Eve 
different cuts on the magnitude of the di-lepton rapidity, \yii\. This time, we consider the 
average over up and down-quarks. From here, one can see that, in searching for extra Z'- 
bosons with masses larger than 0(3 TeV) the \yii\ cut is not mandatory. The true direction 
of the quark is indeed correctly guessed more than 70% of the times, even if no cut is applied 
on the di-lepton rapidity. This means that, at high di-lepton invariant masses, we should be 
able to observe a lepton asymmetry with a well approximated shape even without imposing 
ad hoc cuts. The advantage would be twofold: preserving a small statistical error on that 
shape, owing to the much larger acceptance one should have in absence of the \yij\ cut, 
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FIG. 7: (Colour online) Upper plot: Reconstructed forward-backward asymmetry as a function of 
the di-lepton invariant mass within the SM at the 13 TeV LHC with a total integrated luminosity 
C = 300/6“^ for a set of different rapidity cuts on the di-lepton system. In the legend, \yee\ 
corresponds to \yij\ defined in the text. The black line represents the true AFB for comparison. 
Lower plot: Acceptance as a function of the di-lepton invariant mass for the same set of di-lepton 
rapidity cuts as above. 


and working with an event sample flavour independent up to a large extent. This latter 
feature would guarantee a more model independent procedure, as the different Z' models 
have obviously different couplings of the extra gauge boson to up and down-quarks. 

In order to quantify the delicate balance between AFB line shape and statistical error, 
in the upper plot of Fig. 7 we show the shape of the reconstructed lepton asymmetry, 
within the SM as a function of the di-lepton invariant mass for a set of different cuts on 
\yil\. We compare the results to the true AFB, where the direction of the valence quark is 
taken directly from the Monte Carlo (MC) event generator. In the lower plot of Fig. 7, we 
display the acceptance as a function of the same variable Mu for the same set of \yii\ cuts. 
Comparing the two plots, one can see that A*p^ tends to the true AFB with increasing the 
\yij\ cut, but at the same time the acceptance heavily decreases. In particular, for masses 
above 2.5 TeV, if we apply the stringent cut \yii\ > 0.8 used in literature, the number of 
events goes down by a factor of 3 while the gain in shape is only about 20% of the true 
AFB value. With increasing mass, the acceptance decreases indeed more rapidly with the 
\yii\ cut. 

To visualize how the above features impact the AFB sensitivity to new physics, in Fig. 8 
we compare the reconstructed A*pQ observable predicted by three representative Z'-models 
El and GLR — LR) with the SM expectation at the 13 TeV LHC with total integrated 
luminosity C = 100/6“^. As a new physics signal, we consider a hypothetical Z'-boson with 
mass Mz' = 3 TeV. To quantify the effect of the di-lepton rapidity cut on the signihcance 
either in searching for new physics via AFB or in distinguishing between different Z' models, 
we show results for the commonly used \yij\ > 0.8 setup (Fig. 8(d)) versus the > 0.4 
(Fig. 8(c)) and no cut (Fig. 8(b)) scenarios. We further display, in Fig. 8(a), the ideal 


situation represented by the true forward-backward asymmetry, Afb- 

As one can see, imposing a strong di-lepton rapidity cut helps in recovering the true shape 
and magnitude of the forward-backward asymmetry. However, the consequent decrease 
of the number of events is so substantial that the signihcance diminishes drastically with 
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FIG. 8: (Colour online) (a) True forward-backward asymmetry as a function of the di-lepton 
invariant mass as predicted by the SM (black), the (orange), the Ej (magenta) and the GLR — 
LR (purple) models for a Z'-boson with mass Mz' = 3 TeV. The results are for the LHC at 
^/s = 13 TeV and L = 100/6“^. Lower plot: the significance in distinguishing models is displayed. 
The double colour in each bin visualizes the two compared models. El Reconstructed forward- 
backward asymmetry as a function of the di-lepton invariant mass as predicted by the SM (black), 
the Ey^ (orange), the Ej (magenta) and the GLR — LR (purple) models for a Z'-boson with mass 
Mz' = 3 TeV. The results are for the LHC at ^/s = 13 TeV and £ = 100/6“^. No cut on the 
di-lepton rapidity is imposed: \yii\ > 0. Lower plot: the significance in distinguishing models is 
displayed, (c) Same as plot (b) with \yii\ > 0.4. El Same as plot (b) with \yif\ > 0.8. 


increasing the \yif\ cut. Moreover, the implementation of the \yij\ cut accentuates the flavour 
dependence of the results or, in other words, the model dependence of the analysis. As the 
probability of guessing the correct direction of the quark in the reconstruction procedure of 
the AFB as a function of the \yii\ cut depends on the type of quark (up and down-quarks 
react differently to the cut), the reconstructed AFB shows an increased model dependence 
in its response to the \yif\ cut. To exemplify this concept, let us take the third bin from 
the left in Fig. There (Fig. 8 (a)), the Ey and Ej models are degenerate as far as the 
true asymmetry is considered. When we compare the reconstructed asymmetry, we see that 
the two models are not degenerate any more in that bin. The splitting increases with the 
\yij\ cut, as the two models react differently to such a cut, having different couplings of the 
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corresponding Z'-boson to up- and down-quarks. In order to minimize the presence of model 
dependent elements in the analysis, it is thus advisable not to include the di-lepton rapidity 
cut. Hence, in the following, we will be working in a setup where we do not impose such a 
restriction. 


IV. THE ROLE OF AFB IN Z’ SEARCHES: NARROW HEAVY RESONANCES 

The AFB is the observable where the effects of the interference between new physics and 
SM background are maximal. In the Drell-Yan processes, these effects are of course present 
also in the total cross section. They are readily seen in both cases in the di-lepton invariant 
mass. As mentioned repeatedly, constraining the search window for new physics within the 
interval \Mii — Mz'\ < 0.05 x Zlhc guarantees that hnite width and interference effects are 
below the 0(10%) level when compared to the complete new physics signal. Such effects are 
instead an intrinsic part of the AFB and dominate its dynamics. For such a reason, the AFB 
is an intrinsically model dependent variable and in literature has therefore been traditionally 
considered for disentangling different models predicting a spin-1 heavy neutral particle. Its 
role has therefore been cornered so far to the interpretation of a possible Z'-boson discovery 
obtained via the default bump search. 

In this paper, we aim to show that AFB can also be used for searches, directly, as a 
primary variable alongside the cross section itself. In this section we focus on Z'-bosons 
characterized by a narrow width. This is the most common kind of particle predicted 
by theories with an extra U'{1) gauge group. This is also the scenario mostly studied in 
literature. The experimental searches for such an object are tailored on this expectation 
and the corresponding results coming from the data collected at the 8 TeV LHC have been 
summarized in Section H. With respect to the ‘AFB search’, the Z' models can be divided 
into two categories: Z' models with AFB centred on the Z'-boson mass and Z' models with 
shifted AFB. In the next two subsections, we discuss their properties in turn. 


A. Z’ models with AFB centred on peak 


In this subsection, we discuss models where the AFB is peaked on the Z'-boson mass. 
These models belong to the Eq class of theories which predict new narrow width spin-1 
resonances. In the literature, it is known that such models contain one extra neutral gauge 
boson whose width cannot exceed a few percent of its mass: Tz'/Mz' < 5%. Even the 
inclusion of new Z'-boson decay channels into exotic states would not change this estimate. 

We hrst compare the shape of the AFB distribution as a function of the di-lepton invariant 
mass. Mil, with the differential cross section in the same variable. In Fig. 9, we show results 
for two representative Eq models: E^ and Ej. We consider a hypothetical Z'-boson with 
mass Mz' = 3 TeV. In Figs. 9(b) and |9(d) , we display both the true and the reconstructed 


AFB within the two chosen Eq models with and without taking into account the interference 
between the extra Z'-boson and the SM background. As one can see, the role played by the 
interference is extremely important. The AFB shape is drastically modihed by getting its 
peak heavily accentuated. In contrast, the invariant mass distribution is almost interference 
free if the \Mii — Mz'\ < 0.05 x Zlhc cut is imposed, as shown in Figs. 9(a) and 9(c) for 


the two representative Eq models. In interpreting the experimental data coming from AFB 
measurements it is then mandatory to include the interference, no matter what kinematical 
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FIG. 9; (Colour online) (a) Differential cross section as a function of the di-lepton invariant mass 
as predicted by the model for a Z'-boson with mass Mz' = 3 TeV. The results are for the LHC 
at \/s=13 TeV. (b) Reconstructed forward-backward asymmetry as a function of the di-lepton 
invariant mass as predicted by the model for a Z'-hoson with mass Mz' = 3 TeV. The results 
are for the LHC at \/s=13 TeV. No cut on the di-lepton rapidity is imposed: \yii\ > 0. (c) Same 
as plot (a) within the Ej model. |(d)| Same as plot (b) within the Ej model. 


cut is applied. 

In terms of significance, the search for narrow width Z' models with AFB centred on the 
Z' mass is summarized in Fig. 10 for the two representative models and Ei. Within the 
E^ model, the true AFB would give rise to a significance slightly lower than that one coming 
from the usual bump search, as shown in Fig. 10a. The reconstruction procedure of the AFB 
depletes this result but still the two signihcances from cross section and A*p^ are comparable 
over the full di-lepton invariant mass range . Fig. 10b shows that the Ej model is more 
accessible through the AFB than the cross section. There, indeed, the signihcance from 
the true AFB is a factor two bigger than the significance coming from the bump search. 
Once again, the AFB reconstruction pollutes the ideal result. The signihcance from the 
reconstructed AFB gets reduced, but its value remains anyhow only slightly lower than that 
one coming from the resonance search. The Ej model is not unique in this respect, also the 
Es model shares the same property. 

Similar trends are shown by all models belonging to the Eq class of theories and do 
not change when a more realistic setup is considered. Implementing the acceptance cuts 
extracted by the CMS analysis at the 8 TeV LHC {Pt{^) > 25 GeV and 7 ]{l) < 2.5 with 
I = e,p), the shape of the reconstructed AFB, including error bars would in fact 
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FIG. 10: (Colour online) |(a)| Binned significance as a function of the di-lepton invariant mass as 
predicted by the model for a Z'-boson with mass Mz' = 3 TeV. The red line represents the 
significance corresponding to the invariant mass distribution. The blue and green lines show the 
significance extracted by an ideal measurement of true and reconstructed AFB, respectively. The 
results are for the LHC at y/s = 13 TeV and C = 300/6“^. 


(b) Same as (a) for the Ej model. 


appear as in Fig. 11. The significances coming from the AFB and the cross section are 
indeed equivalent in magnitude, if only the statistical error is included. We thus expect 
that the use of the A*p^ observable, when associated to the default resonance search, could 
improve the discovery potential of new narrow width Z'-bosons. Further, being a ratio of 
differential cross sections, the reconstructed A*p^ could help in minimizing the systematical 
errors thus rendering the measurement much more accurate. 

This is in particular the case when we should be in presence of an evidence for a new 
Z'-boson in the resonance search at the 3-4 sigma level. In these conditions, one couldn’t 
claim the discovery of a new gauge boson just looking at the resonant peak in the di-lepton 
invariant mass distribution. However, if a signal of similar strength were to be discovered 
in an independent observable, the suggestion of the possible presence of new physics would 
turn into a robust claim. This is the role that the AFB would play. In Fig. 12, we plot the 
differential cross section and A*pp as a function of in the di-lepton invariant mass, = a/I, 
within the and Ej models, at the forthcoming Run II of the LHC at 13 TeV with 
C = 30 fb~^, that is the integrated luminosity which should be collected at the end of 2016. 
There, a new physics evidence at barely 4 sigma in the bump search could be reinforced by 
the simultaneous measurement of the reconstructed AFB, showing a signal at the 2-sigma 
level. 


B. Z’ models with shifted AFB 

In this section, we discuss narrow width Z' models where the AFB has a shifted peak, 
that is, not centred on the Z'-boson mass. These models belong to the GLR class. The 
same behaviour is also displayed by the SSM scenario taken as benchmark model by the 
LHC experimental collaborations. 

In principle, the reconstructed A*pp could reveal the presence of a new spin-1 particle at 
energy scales lower than its mass, as the shape of this observable as a function of the di-lepton 
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FIG. 11: (Colour online) [(a)] Binned differential cross section as a function of the di-lepton invariant 
mass within the model for a Z'-boson with mass Mz' = 3 TeV. Error bars are included. The 
results are for the LHC at \/s=13 TeV and C = 300/6“^. Acceptance cuts are imposed (see text). 
The lower plot shows the signihcance. |(b)| Binned A*p^ as a function of the di-lepton invariant 
mass within the Ey- model for a Z'-boson with mass Mz' = 3 TeV. Error bars are included. The 
results are for the LHC at y/s = 13 TeV and L = 300/6“^. Acceptance cuts are imposed (see text). 
In the lower plot, the blue histogram shows the binned significance while the green area indicates 
the total significance integrated over that invariant mass region, (c) Same as (a) for the Ej model. 
|(d)| Same as (b) for the Ej model. 


invariant mass is accentuated at mass scales smaller than Mz'- This behaviour is shown in 
Fig. 13b where we plot the reconstructed AFB versus Mii = \/J for the representative model 
GLR-LR. We consider a new Z'-boson with mass Mz' = 3 TeV. As one can see, the peak of 
A*pp is shifted on the left-hand side of the physical Z'-boson mass and it appears at around 
2.6 TeV. This feature is interesting. However, the signihcance is quite low as shown in Fig. 
13b, owing to the poor statistics in that region. For Mij values around the physical mass of 
the Z'-boson, which are statistically relevant, the signihcance coming from AFB is always 
much smaller than the signihcance obtained via the measurement of the differential cross 
section, displayed in Fig. 13a. For this kind of models, the A*pp observable is therefore not 
particularly appropriate for Z' searches. The same conclusion holds for the SSM (see Figs 
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FIG. 12: (Colour online) [(a)] Binned differential cross section as a function of the di-lepton invariant 
mass as predicted by the Ey, model for a Z'-boson with mass Mz' = 3 TeV. Error bars are included. 
The results are for the LHC at \/s=13 TeV and C = 30/5“^. Acceptance cuts are imposed (see 
text). The lower plot shows the significance. |(b)| Binned A*p^ as a function of the di-lepton 
invariant mass as predicted by the Ey. model for a Z'-boson with mass Mz' = 3 TeV. Error bars 
are included. The results are for the LHC at \/s=13 TeV and L = 30/6“^. Acceptance cuts are 
imposed (see text). In the lower plot, the blue histogram shows the significance bin by bin while 
the green area indicate the total signihcance integrated over that invariant mass region, (c) Same 
as (a) for the Ei model. |(d)| Same as (b) for the Ei model. 


13c and 13d). Hence, this benchmark model is not an advisable playgronnd for stndying the 
benehts of nsing the AFB in searching for new Z'-bosons. 


V. THE ROLE OF AFB IN Z’ SEARCHES: WIDE HEAVY RESONANCES 

In this section, we discnss the role of the reconstrncted AFB, in searches for a new 
Z'-boson characterized by a large width. Snch a heavy and wide particle is predicted by 
different models. A benchmark scenario for experimental analyses is the wide version of the 
SSM described in Ref. [Sj. The proposal is to have a heavy copy of the SM nentral gange 
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FIG. 13: (Colour online) (a) Binned differential cross section as a function of the di-lepton invariant 
mass as predicted by the GLR — LR model for a Z'-boson with mass Mz' = 3 TeV. Error bars 
are included. The results are for the LHC at t/s= 13 TeV and L = 300/6“^. Acceptance cuts are 
included (see text), (b) Binned A*p^ as a function of the di-lepton invariant mass as predicted by 


the GLR — LR model for a Z'-boson with mass Mz' = 3 TeV. Error bars are included. The results 
are for the LHC at -v/s=13 TeV and £ = 300/5“^. Acceptance cuts are included (see text), (c) 
Same as (a) for the Generalized Sequential SM, GSM — SM. |(d)| Same as (b) for the GSM — SM 
model. 


boson, Z, with same couplings to ordinary matter and SM gauge bosons. Owing to the 
Z'-boson decay into SM charged gauge bosons, whose rate grows with the third power of the 
Z'-boson mass, the total width of the new heavy particle can be quite large: Tz'/Mz' — 50% 
and above. 

In this case, the invariant mass distribution of the two hnal state leptons does not show 
in the cross section a resonant (or peaking) structure around the physical mass of the Z'- 
boson standing sharply over a smooth background, but just a shoulder spread over the SM 
background. This result is plotted in Fig. 14a, where we consider a Z'-boson with mass Mz' 
= 1.5 TeV and width Vz' jMz' = 80%. As the line shape of the resonance is not well defined 
and these parton level results could be worsened by detector smearing effects giving rise to 
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FIG. 14: (Colour online) (a) Binned differential cross section as a function of the di-lepton invariant 
mass as predicted by the GSM — SM model for a Z'-boson with mass Mz' = 1.5 TeV and 
Vz'lMz' = 80%. Error bars are included. The results are for the LHC at \/s=13 TeV and 
£ = 300 / 6-1 


Acceptance cuts are included (see text), (b) 


Binned A*p^ as a function of the 
di-lepton invariant mass as predicted by the GSM — SM model for a Z'-boson with mass Mz' = 
1.5 TeV and Vz' jMz' = 80%. Error bars are included. The results are for the LHC at y^=13 TeV 
and £ = 300/6-1. Acceptance cuts are imposed (see text). 


an even broader spectrnm, the A*pp observable conld help to interpret a possible excess of 
events. The resnlts are shown in Fig. 14b. There, one can see that the A*pp shape could be 
visible at the 2cr level. 

A framework, theoretically more grounded than the wide SSM, which predicts a heavy 
and broad Z'-boson is the so-called non-universal SU{2) model [SI Ej- In this theory, the 
third generation of fermions is subjected to a new SU{2) dynamics different from the usual 
weak interaction advocated by the SM. On the contrary, the first two families of fermions 
only feel the SM weak interaction. As a consequence, a new spectrum of gauge bosons 
emerges in the model. These new vector bosons can be either narrow or wide. The only 
constraint comes from the EW Precision Tests (EWPTs) which bound the Z'-boson to have 
a mass Mz' > 1.7 TeV. For the analysis of the direct and indirect limits on this model we 
refer to and references therein. 

Within this framework, we consider the wide Z'-boson case with Yz' jMz' = 50%. In order 
to fulhl both the limits quoted in Ref. [7| and the direct limits coming from direct searches 
at the 8 TeV LHC [13], we assume Mz' = 5.5 TeV. The latter analysis performed at the 
LHC has been optimized for searches of new physics with no resonant peaking structure. 
The outcome is that there are no events for di-lepton invariant masses above 1.8 TeV. We 
have taken this limit into account when evaluating the Z'-boson mass and width. 

This model is a very good playground to test whether the AFB can be used as a primary 
variable in searches for wide objects. In this case, in fact, the new physics signal appears 
as an excess of events spread over the SM background. Almost no line shape is present in 
the di-lepton invariant mass distribution usually measured. Searches are performed relying 
on a pure counting strategy, a procedure which does not allow much interpretation of the 
hypothetical signal. The exploitation of the reconstructed A*pp could help in this respect. 
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FIG. 15: (Colour online) (a) Binned differential cross section as a function of the di-lepton invariant 
mass as predicted by the TH model for a Z'-boson with mass Mz' = 5.5 TeV and Tz' jMz' — 50%. 
Error bars are included. The results are for the LHC at \/s=13 TeV and C, = 300/6“^. Acceptance 
cuts are imposed (see text), (b) Binned A*p^ as a function of the di-lepton invariant mass as 
predicted by the TH model for a Z'-boson with mass Mz' =5.5 TeV and Tz' jMz' — 50%. Error 
bars are included. The results are for the LHC at y^=13 TeV and L = 300/6“^. Acceptance cuts 
are imposed (see text). 


In Fig. 15, we compare the Z'-boson spectrum (15a) and the reconstructed A*pp (15b) 
as functions of the di-lepton invariant mass. As one can see in Fig. 15a, the cross section 
spectrum at parton level is already very broad. Its slope might be lost or mistaken in the 
SM background normalization. Even if, in the best case, a plateau would be visible over 
the SM background, its interpretation would be very difficult. Fig. 15b shows that the 
A*pb observable has a sharper line-shape which can reveal the presence of a spin -1 particle 
beyond error bars. Such a shape is quite shifted at low energy scales though, compared 
to the Z'-boson mass. Hence, its extraction should enable one to help the discovery of a 
new vector boson with very high mass. In short, here, the AFB measurement could become 
particularly useful at the edge of the LHC discovery limits, when new particles can be too 
heavy and broad to be easily detected via a standard resonant peak search. 

The aforementioned scenarios are in fact particularly challenging for experimentalists. 
The non-resonant analyses of wide objects have been performed by searching for a smooth 
deviation from the SM background. The number of events above a given lower cut on the 
di-lepton invariant mass is compared with the total number of expected background events. 
An optimal minimum mass threshold is chosen to maximize the sensitivity to new physics. 
Clearly, such an analysis depends quite strongly on the SM background estimate. Usually, 
the simulated background is normalized to the event number in a mass window of ± 30 GeV 
around the Z-boson mass. A control region is then selected at higher di-lepton invariant 
masses in order to perform a data driven modelling of the SM background and recast it 
in a functional form easy to implement in the likelihood used for extracting the limit on 
the Z'-boson mass. The method is based on the assumption that the control region is new 
physics free. But, this is not the case for wide Z'-bosons. In these scenarios, the interference 
between the extra Z'-boson and the SM 7 , Z is so sizeable that it can invade the control 
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region. Being absolutely model-dependent, it can be either constructive or destructive. In 
any case, it would change accordingly the shape of the di-lepton spectrum. If the interference 
is negative, it would led to a depletion of events at low mass scales on the left-hand side of 
the Z'-boson resonance. This is exactly the example shown in Figs. 14 and 15 corresponding 
to the SSM and non-universal SU{2) scenarios, respectively. If not correctly interpreted, 
these interference effects could induce one to underestimate the SM background with the 
consequence of overestimating the extracted mass bounds. Having all these uncertainties to 
deal with, the support of a second observable like the AFB is strongly recommended for the 
non-resonant analyses. 


VI. ON THE ROBUSTNESS OF AFB AGAINST PDF UNCERTAINTIES 

In this section, we discuss how robust is the shape of the forward-backward asymmetry 
against the theoretical uncertainties on the PDFs. We further compare the systematic error 
induced by the PDFs uncertainty on the differential cross section and on the reconstructed 
AFB. 

For the determination of the PDF uncertainty we follow [22] and references therein. 
Here, we just highlight the key points of the procedure. We compute the Hessian PDF 
uncertainty for our two observables: di-lepton invariant mass distribution in cross section 
and reconstructed AFB. For Hessian PDF sets, both a central set and error sets are given. 
The number of error sets is twice the number of eigenvectors. For the CTEQ6.6 PDF that 
we use, the number of error sets is equal to 40. Dehning Xf the value of the variable using 
the PDF error set corresponding to the ”±” direction for the eigenvector i, the symmetric 
error on the variable X is given by: 
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With this dehnition, we are now ready to compute the PDF uncertainty of any function 
f{X). For the differential cross section, we apply Eq. [^directly. For the AFB, the com¬ 
putation is slightly more involved. We consider as independent variables the forward and 
backward (differential) cross sections: ap and ap- The PDF degrees of freedom of these two 
observables are correlated, that is the quantity 
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is equal to 1. Under this condition, by applying the error chain rule, we get 


AApB ~ 


.AcTp Acjf 
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( 10 ) 


( 11 ) 


The sign appearing in the above formula is crucial for the AFB. It indeed clearly shows 
that there is a partial cancellation of the PDF error on the reconstructed App due to the fact 
that this observable is a ratio of (differential) cross sections. Compared to the differential 
cross section, the AFB is then more robust against PDF uncertainties. This is shown in Fig. 
16 where we compare the effect of PDF and statistical errors on the shape of the di-lepton 
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FIG. 16: (a) Differential cross section as a function of the di-lepton invariant mass as predicted 


by the model for a Z'-boson with mass Mz' = 3 TeV. The results are for the LHC at \/s=13 
TeV and L = 300/6“^. The solid line shows the central value, the dotted line the PDF uncertainty. 
The inset plot displays the ratio between PDF and statistical errors, (b) A*p^ as a function of the 


di-lepton invariant mass as predicted by the model for a Z'-boson with mass Mz' = 3 TeV. The 
results are for the LHC at -y/s=13 TeV and C = 300/6“^. The dotted lines show the PDF error 
band, while the dashed lines define the statistical error band, (c) Same as (a) for the Ej model. 
|(d)| Same as (b) for the Ej model. 


invariant mass distribution of the cross section and AFB for two reference models, Ej and 

As one can see, the behaviours of cross section and AFB are opposite. The differential 
cross section in the di-lepton invariant mass is dominated by the PDF error on and off-peak. 
In the region around the peak, the PDF uncertainty is a factor 2 bigger than the statistical 
error. On the contrary, the AFB is dominated by the statistical error on and off-peak. 
Moreover, the PDF uncertainty is quite reduced owing to the minus sign in EqjTTj The 
shape of the AFB is thus not affected by the PDF error, so this observable is theoretically 
well dehned. 

In the light of these results, we can revisit Figs 11 and 12. We immediately observe 
that, if the error from the PDF is a factor 2 bigger than the statistical one for the invariant 
mass distribution around the Z'-boson mass, then the total error becomes 3 times the 
statistical one (adding linearly the systematic and statistical errors). As a consequence. 
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FIG. 17: (Colour online) (a) Cross section integrated around the Z'-boson mass — Mz'\ < 
0.05 X F'lhc) as a function of Mz' as predicted by the E^, Ej, GLR — LR and GSM-SM models. 
The results are for the LHC at -v/s=13 TeV and £ = 300/6“^. Solid lines represent the PDF 
uncertainty, dashed ones the statistical error, (b) A*p^ integrated around the Z'-boson mass 
{\Mij — Mz'\ < 0.05 X F'lhc) as a function of Mz' as predicted by the E^, Ej, GLR — LR and 
GSM-SM models. The results are for the LHC at v^=13 TeV and C = 300fb~^. Solid lines 
represent the PDF uncertainty, dashed ones the statistical error. 


the significances in Figs 11a and 11c decrease by the same factor. This means that the 
signihcances from the corresponding AFB distribntion are mnch bigger than those ones 
from the differential cross section. If no rehtting procednre is applied for the PDF in order 
to minimize their nncertainty, then the AFB observable appears to be mnch more robust 
than the resonance peak. 

Of course, one needs to consider the energy scale dependence of the PDF errors, if no 
rehtting procedure is employed. In Fig. 17a, we plot PDF and statistical errors on the 
total cross section integrated around the mass of the Z'-boson. We integrate in the window 
±5% F'lhc around the hypothetical Mz> where interference and hnite width effects can be 
neglected. In evaluating the statistical error, we assume the design value for the luminosity: 
L = 300/6“^. We then vary the value of Mz' to see how statistical and PDF errors change 
in magnitude. We consider four theoretical frameworks: E^, Ej, GLR-LR and SSM. The 
hgure shows that, up to roughly a 4 TeV scale, the cross section is dominated by the PDF 
uncertainty. In contrast, the asymmetry integrated in the same peak region is heavily 
dominated by the statistics for all possible Z' masses, as shown in Fig. 17b. 

The strong dependence of the PDF errors on the mass or energy scale also suggests that 
using as observable the ratio between the Z'-boson cross section and the on-peak SM Z- 
boson cross section R„, might not be entirely PDF safe. The two cross section are indeed 
a few TeV a part. Given the strong variation of the PDF error with energy, a cancellation 
analogous to that one in Eq. 11 could not happen easily. This is another argument in favour 
of exploring the AFB as a search variable. 


VII. CONCLUSIONS 

In this paper we have considered the scope of using AFB, the forward-backward asymme¬ 
try, in Z'-boson searches at the LHC in the di-lepton channel, i.e., via Drell-Yan production 
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and decay. Such a variable has traditionally been used for diagnostic purposes in presence of 
a potential signal previously established through a standard resonance search via the cross 
section. However, based on the observation that it is affected by systematics less than cross 
sections (being a ratio of the latter), we have studied the possibility of using AFB for such a 
purpose for a variety of Z' models, GLR, GSM, embedding either a narrow or wide res¬ 
onance. The focus was on determining whether such a resonance could be sufficiently wide 
and/or weakly coupled such that a normal resonance search may not fully identify it and, 
further, whether the AFB could then provide a signal of comparable or higher signihcance 
to complement or even surpass the scope of more traditional analyses. 

We have found promising results. In the case of narrow width Z'-bosons, we have proven 
that the signihcance of the AFB search can be comparable with the usual bump search. 
Further, we have emphasised the fact that the AFB distribution mapped in di-lepton in¬ 
variant mass can present features amenable to experimental investigation not only in the 
peak region but also signihcantly away from the latter. In the case of wide Z'-boson, the 
AFB search could have a better sensitivity than the cross section studies thanks to a more 
peculiar line-shape and lower systematic and PDF uncertainties. In essence, here, AFB in 
specihc regions of the invariant mass of the reconstructed Z'-boson could be sensitive to 
broad resonances much more than the cross section, wherein the broad distribution of the 
signal seemingly merges with the background. 

We have explored the above phenomenology for all the benchmarks under study as well 
as assessed and used the optimised strategy for AFB reconstruction. 
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